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INTRODUCTION 

A limited number of small, heavy-metal-reflected, fast reactor, 
critical experiments have been performed at Los Alamos (Reference 8). 
These have been used to check the accuracy of neutronic calculations. 
From these calculations it has been found that the amount of c&lcula- 
tional error is strongly dependent on the materials in the reactor and 
in the neutron reflector. 

Since for the materials of interest the available experiments 
supply only about one data point, and there are some materials such as 
nitrogen and hafnium for which there are no fast spectrum critical ex- 
periments, it was decided to embark on a critical experiment program 
wherein only the materials that were being considered for an actual, 
small, fast- spectrum reactor would be included. Also, in order to 
minimize the number of geometrical variables while determining the re- 
activity effects of the proposed materials, it was decided to make an 
initial critical experiment large enough and with sufficient void 
space so that all of the materials could be added sequentially without 
changing the overall geometry of the reactor. This critical experi- 
ment program is being carried out by Atomics International under con- 
tract to LeRC. (Contract No. NAS3-12982) . 

The conceptual reactor for which these experiments are being per- 
formed is a lithium- cooled , uranium-235 nitride fueled reactor. A cut- 
away drawing is shown in Fig. 1. Cylindrical uranium nitride pellets 
would be contained in tantalum alloy (T-lll) tubes which would be lined 
with tungsten foil. A T-lll honeycomb would be used to maintain the 
position of the fuel elements. Six control drv r s containing both fuel 
and neutron absorber (T-lll) would be used for reactivity control. The 
neutron reflector would be of a molybdenum alloy (TZM). 

Since there are many ways (voluntary and involuntary) in which a 
neutronic calculation can be made to agree with a critical experiment 
after the experimental results are known, and since for a final reactor 
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design one would like to be able to truly predict the core parameters, 
the purpose of this report is to present a precritical analysis of the 
fuel loading requirements and neutron energy spectrum of the first of 
a series of critical experiments which are to be performed by Atonies 
International. 


Description of the Critical Assembly 

A cross section of a fuel element is shown in figure 2. It is 
depicted as filled with a full complement of fuel which is oralloy 
(93*2 percent uranium-235 with the remainder assumed to be uranium-238). 
The space between the fuel bundle and the fuel tube can be used to insert 
foils of other materials such as tungsten, hafnium, and tantalum. Both 
the fuel tube and honeycomb tube are made of tantalum. The fuel tube 
is centered offset from the center of the honeycomb tube by 0.020*0.002 
inch (0.051*0.005 cm) so that a fuel movement experiment can be performed 
simply by rotating the fuel elements. The eccentric annular space be- 
tween the fuel tube and honeycomb tube can also be used for addition of 
materials of interest such as lithium nitride to simulate the lithium 
which, as lithium 7, may be used for coolant in a power reactor and 
nitrogen to simulate the ni'rogen in uranium nitride fuel. 

The end reflectors are built into the fuel elements as shown in 
figure 3> which shows one end of a typical fuel element. In each fuel 
element there are two molybdenum pieces at each end of the fuel element. 
One piece is a solid cylinder 3. 94 inches (10.00 cm) long which fits 
within the fuel tube and rests against the ends of the 14. 767 -inch 
(37* 508- cm) long fuel bundle. The other piece fits in the eccentric 
annular space between the fuel tube and honeycomb tube, bume support 
must be provided for this piece to prevent its sliding into the core 
region. Several small diameter molybdenum rods the same length as the 
fuel bundle could be used. For the purpose of the calculations presented 
later, three 0.030-inch (0.076-cm) diameter rods per fuel element are 
used. 


The ends of the fuel element are closed with an aluminum Jacks crew 
and rubber packing as shown in figure 3* The grid plate alinement pin 
provides a method of constructing a lattice of these fuel elements with 
a constant pitch by the use of top and bottom grid plates. Figure 4 
shows a cross section of a quadrant of the critical assembly. Only the 
outline of the honeycomb tubes are shown to depict the lattice deter- 
mined by the grid plate. The center-to-center spacing (pitch) between 
tubes in the nondrum region of the core is O .872 inch (2.215 cm). 
Similarly grid plates for each control drum position the fuel elements 
as shown. The assembly consists of 24? fuel elements including 11 in 
each control drum. The individual pieces called out on figure 4 are 
described in more detail in figures 5 and 6 in table I. Figure 5(a) 
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shows the molybdenum filler pieces that fit between the rotatable con- 
trol drums and the core lattice. Figure 5(b) shows the radial re- 
flector pieces that fit between the control drums. In the calculations, 
these pieces are assumed to be pressed and sintered with a density of 
0.9 h of theoretical (10.22 gm/cnr). Four pieces as shown are used. 

Two of the pieces are turned upside down so that the arrangement of a 
pair of reflector pieces about a control drum as shown in figure 4 can 
be obtained. Two additional reflector pieces are required. These are 
the same as the others except that the piece is symmetric with the 
shape of side A (fig. 5(b)). Figure 6 shows a control drum in more de- 
tail. Table I gives additional data on individual pieces used in the 
reactor. Figure 7 shows a sketch of the assembled critical machine. 

This sketch was supplied by the contractor. Note that the safety scram 
device consists of moving two pairs of reflector pieces away from the 
core and that the machine is elevated to provide for a sample changer 
for the insertion of materials for reactivity measurements. The sam- 
ples may be inserted into the center seven fuel element locations or 
into the six positions outer-most in the core lattice. 


Calculation Procedure 

The calculations for this critical assembly uses basically the 
same techniques as described in references 1 to 3* As stated previously 
the purpose of these precritical calculations is to provide evidence as 
to the adequacy of these methods. Discrepancies, if any, between the 
experimental quantities and the calculations will be due to inadequate 
methods, inadequate cross section data, differences between the calcula- 
tional geometric model and the "as-built" critical assembly, or any 
combination of these. A detailed post-critical analysis of the "as- 
built" critical assembly is desirable to evaluate all of the parameters 
used in the calculatlonal procedure. 

The calculations reported here are based upon the nominal dimen- 
sions and material specifications in the latest drawings available at 
this time. Three types of spatical calculations are used to obtain the 
expected critical mass of oralioy fuel. A two-dimensional (X,Y) calcu- 
lation with the mesh grid superimposed over a quadrant of the reactor 
as shown in figure 4 allows explicit representation of the azimuthal 
asymmetry around the periphery of the core. An alternate calculation 
that is also suitable uses (R,9) geometry; both (X,Y) and (R,©) calcu- 
lations are done for this critical experiment geometry. The geometry 
in the direction perpendicular to the (X,Y) or (R,©) plane is accounted 
for by inclusion of an effective core height He. f from which axial 
neutron leakage terms may be computed. The procedure for determining 
HOff is as follows: 
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1. The (X,Y) or (R,0) geometry is converted into a series of 
annuli by volume averaging materials in the two-dimensional geometry 
that appear in these annuli. An (R,Z) model, such as in figure 7, is 
established which shows the core height and axial reflector regions 
explicitly within the restrictions implicit in the (R,Z) synmetry re- 
quirement. 

2. A one-dimensional (R) calculation, based on the radial di- 
mensions and materials of the (R,Z) model, is used to obtain He-~ • He- 
in the formula for the axial leakage probability per unit flux 



is varied until the multiplication factor from the R calculation is 
the same as from the (R,Z) calculation. This He— is then used in 
the (R,©) or (X,Y) calculation. In the above formula is the 

transport cross section for group g and material i // . . 
in the R calculation; B 2 ° ir/y37 g 

The GAM-II program (Ref. 4) is used to obtain the cross sections. 
Three spectra are used in the averaging of the bro*vi group microscopic 
cross sections. In each of the three cases, the spectrum is obtained 
allowing a uranium-235 fission spectrum to interact with the homogeneous 
materials. Both 13 group and 4 group cross sections with the energy 
groups as shown in table H are obtained. The spectra used in the aver- 
aging are representative of the homogeneous core, molybdenum radial re- 
flector and homogeneous axial reflector regions. Additional materials 
appearing in these regions, or nearby, vnich would not contribute 
significantly to the spectrum are nevertheless averaged over the spectrum 
for use in the spatial calculations, as an example, the pressure vessel 
tantalum is averaged over the radial molybdenum reflector spectrum. The 
MACROS program is used to produce macroscopic cross sections for use in 
TD8N (Ref. 5 ) or DOT (Ref. 6) Which are used for spatial solutions. 


RKSULT8 

Results of the analysis consist of the spectrum of neutrons at the 
center of the core, the prediction of critical mass, and some additional 
reactivity values calculated to determine the sensitivity of the multi- 
plication factor to perturbations to the basic system either physically 
or strictly from a computational viewpoint. 
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Flux 8pectrum 

The flux spectrum Is obtained from the GAM-II computer program in 
order to show the spectrum in more detail than is possible with the 
broad energy group spatial solution programs. A homogeneous mixture 
of fuel cell materials (0.0692 volume percent v/o Ta, 0.2413 v/o or alloy, 
0.0032 v/o Mo) corresponding to a reactor fuel loading of 179*49 kg of 
oralloy was used. Figure 9 shows the resulting spectrum. The median 
energy of this spectrum is O.78 MeV. 


Criticality Calculations 

Calculations for the critical assembly using computer programs and 
techniques currently available show considerable variation in the com- 
puted multiplication factors. This fact argues stx'ongly in favor of a 
definitive critical experiment program to reduce computed uncertainties 
to acceptable values for the design of a nuclear space-power reactor. 
Table HI shows computed multiplication factor ke. _ for a fuel loading 
of 179*49 kg of oralloy. The two-dimensional RZ calculation was per- 
formed to incorporate the effects of the end reflectors and structure 
and to be used in obtaining an effective core height for use in the other 
two-dimensional calculations (R9 and XY) . 

For the XY calculations a mesh grid of 39 intervals in the X 
direction and 35 intervals in the Y direction were used. The pressure 
vessel was not included in these calculations although its worth 
(O.OO56 dk) has been added to the multiplication factors in table III* 

The worth was calculated using the RZ model. With the XY geometry, in- 
creasing the order of the angular quadrature from S« to 8v increases the 
multiplication factor from 1.02 56 to 1.0357* This Increase has been ob- 
served previously for similar reactors though the effect was not as pro- 
nounced. The effect seems related to the mean free path of neutrons in 
the core) this core has long mean free paths due to the largo amount of 
void (about 69 volume percent). Similarly with calculations involving 
curved surfaces such as RZ and R9 gecmetrices the tendency is for the 
multiplication factor to decrease in going from 0« to Sv . Again, the 
magnitude of the decrease is related to the mean Tree path in the core. 
This decrease is seen in table III for the R9 calculations and amounts 
to 0.0059 dk. It is also disconcerting to note the tendency of the R9 
calculation to yield larger multiplication factor than the XY calcula- 
tions. This tendency has also been noted in previous calculations for 
similar reactors. 

The %Po 4 group AY calculations are used to estimate the critical 
loading of the critical assembly as shown in figure 10. The upper curve 
represents the computed values* The lower curve results from sub- 
tracting 3*6 percent dk/k bias from the upper curve. The estimated 
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critical mass of oralloy fuel is thus 179*7 kg. The -3*6 percent 
Ak/k is based on the following: 

1. The use of 4 energy groups instead of 13 or more accounts 
for about -0.4 percent Ak/k. 

2. The use of S^Fo instead of a higher order calculations such 
as SqP^ accounts for about -0.9 percenc Ak/k, 

3. Structure around and at the ends of the reactor, such as drum 
drive motors, scram devices and sample changer mechanist' (see Fig. 6), 
can add about +0,6 percent Ak/k. 

4. Previous calculations of compact critical assemblies fueled 
with oralloy and containing or reflected by refractory metals have, 
in the past, tended to yield multipli cation factors about 2 percent 
Ak/k too large (Ref. 7). Table IV summarizes more recent calculations 
done by R. M. Westfall at this laboratory. These are to be compared 
with an experimental k gff = 1.000. From these it is concluded that the 
undefined part of the calculational bias should be -2.9 percent. The 
four items result in a net 3*6 percent Ak/k that is removed from the 
2HJCY S^Po 4 group calculated values. 

Additional calculations have been made to determine the sensitivity 
of the multiplication factor to perturbations to the critical assembly. 

The results are described in table V. Cases 1 and 2 pertain to the effect 
of additional material around the reactor. Case 3 was calculated be- 
cause the fuel rods are coated with a thin layer of protective material 
containing about 0.4 percent by weight of hydrogen. The atom density 
of hydrogen in a fuel cell is about lxlO 1 ^ atoms/cmS. Case 5 was com- 
puted with the RZ model to use as a correction for the XT model results. 
The relatively large worth of the pressure vessel indicated that materials 
outside the reactor as calculated could be Important (Cases 1 and 2). 

Case 6 related to a way of supporting the upper eccentric molybdenum re- 
flector piece by use of 3 molybdenum rods per fuel assembly between the 
honeycomb tube and fuel tube. Cases 3> 6 and 7 indicates the reactor 
is not overly sensitive to small amounts of material uncertainties in 
the core. Case 8 represents a fuel coefficient in the vicinity of 177 kg 
of oralloy. Cases 9 through 12 were computed using a one-dimensional 
radial calculation and were used in estimating the bias that might be 
expected for the S. Po 4 group XT model used to estimate the critical 
loading. 
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CONCLUDING REMARKS 

An analysis of the fuel loading requirements and. the sutron 
energy spectrum of a critical experiment which is to be per -armed by 
Atomics International under Contract No. NAS3-12982 is presented. 

The critical loading is predicted to be about 179 kg of oralloy. 
The median energy of the neutron energy spectrum is predicted to be 

0.78 MeV. 
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TABLE I. - DESCRIPTION OF VARIOUS REACTOR COMPONENTS 
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TABLE II. - GAM ENERGY GROUP BOUNDARIES 



TABLE III. - CALCULATED MULTIPLICATION FACTORS 


[Fuel loading 179*49 kg cralloy] 
[4 Energy grc p Po calculations] 


Calculation geometry 

program 

Sn 

Multiplication factor 

RZ 

TDSN 

S 4 

l.v.379 

XY 

TDSN 

S 2 

1.0256 

XY 

TDSN 

S 4 

-■.0151 

R0 

DOT 

S 2 

1.0490 

Re 

1 

DOT 

S 4 

1.0431 
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TABLE IV. - CALCULATED MULTIPLICATION FACTORS (k eff ) 


OF CRITICAL EXPERIMENTS 


' — 

Core composition 

Reflector 

Calculation 

keff 

Reference 

50.8 v/o Ta~ Or alloy 

Hone 

S lf P 1 13 group 

1 .C 29 

8 (P. 17) 



1 D synthesis 



Oralloy 

1" Molybdenum 

Vl 13 group 





2 DBZ 

1.029 

8 (P. 17) 

OrrJLloy 

1" Tungsten 

W 13 group 





2 DBZ 

1.007 

8 (P. 28) 


TABLE V. - EFFECT OF H9RTURBATI0HS TO REACTOR OH MULTIPLICATION FACTOR 


- 

r 11 ■ ■— - - -■ ■ — 

Perturbation 

Change in multiplication 
factor, percent 

fl 


40.34 

2 

Add 6-30° sectors of 1.27 cm thick 0.9 
dense stainless steel to pressure vessel 
in area between drums 

+ .31 

3 

Worth of hydrogen in coating on fuel rods 

4 .09 

4 

Core length increased 1 mm 

4 .048 

5 

Removal of pressure vessel 

- .56 

6 

Removal of 3-0-030" dia. Mo rods per 
fuel tube in all fuel tubes 

- .074 

7 

Use of +10 percent tolerance (instead of 
nominal) on fuel and honeycomb 
tantalum tube radial tolerances; 
amount to ~0.7 v/o of a fuel cell 

4 .Olfl 

8 

Add 1 kg of oralloy fuel uniformly 

+ .33 

B 

S^Po 4 group — ^ 8 ^ Po 13 group 


9 

Bo 13 group — b 8g Po 13 group 

- .39 

9 

Sg Po 13 group — > Sg Pj^ 13 group 

- .56 

tfj 

Sg Po 13 group 8^ Po 13 group 

- 09 



















CD-1041?*# 





* n 







Fve\ FUr*** f r r _ 

Zenecroz. o7s 

All dwensionj in inches^ 



//.2£Z 




-7.iUt.fl 


i.ZS 



t 

,38 0 


A* 


FJUt ?‘*ce 


i 


24J+R 

I 


tQ 


if* 

§ 


s 

+t 

* 


fc 


St 
♦ • 


'U/ 



\ 


\ 




f iq.S 

Mx DlW^loWi Ihl INCHES 


~ \ „ 



• !>»*• 


If cf* *yCc/*b *WW& 


° # 

® h ^ 

O 

/ co*4+*! Vtttwi Gtecwel'r/ 

>3 •• i 

Dtta%£*Si**9 /*■? i*ch*>. 















( #cu3t£/ Niuuondvaotf 


M ON XllMt 


NASA-L«wi». 

UVd -XH3 


